Gravitino dark matter and neutrino masses with bilinear R-parity violation 
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Bilinear R-parity violation provides an attractive origin for neutrino masses and mixings. In such 
schemes the gravitino is a viable decaying dark matter particle whose R-parity violating decays lead 
^ , to monochromatic photons with rates accessible to astrophysical observations. We determine the 

^ ■ parameter region allowed by gamma-ray line searches, dark matter relic abundance and neutrino 

■ oscillation data, obtaining a limit on the gravitino mass < 1-10 GeV corresponding to a rela- 
tively low reheat temperature Tr < few xlO'' — 10* GeV. Neutrino mass and mixing parameters 
may be reconstructed at accelerator experiments like the Large Hadron Collider. 

a; 

D . I. INTRODUCTION 

] The origin of neutrino masses and mixing and the nature of dark matter are two of the most elusive open problems 
^ ■ of modern particle physics and cosmology, which clearly indicate the need for new physics beyond the Standard Model. 

■ It has been suggested that these two apparently unrelated issues may be closely inter-linked 

m. 

Here we propose 

■ an alternative way to relate dark matter with neutrino properties within a scenario where supersymmetry (SUSY) is 
. the origin of neutrino mass Q , thanks to the spontaneous violation of R-parity Q . For definiteness and simplicity we 
, adopt an effective description in terms of explicit bilinear R-parity violating superpotential terms (BRpV) [7|-[9| . 

We show how both dark matter and neutrino oscillations can be simultaneously explained in the presence of bilinear 
R-parity violation with gravitino lightest supersymmetric particle (LSP), in such a way that. 



, • gravitino dark matter properties are closely related to the scale of neutrino mass, and 

X' . . . 

5^ ■ • neutrino oscillation parameters may be reconstructed at accelerator experiments. 



Indeed, in this model the very same lepton number violating superpotential terms that generate neutrino masses and 
mixing also induce dark matter gravitino decays, as this also breaks R parity. 

We show how, although unprotected by R-parity, the gravitino can be stable over cosmological times and be a viable 
cold Dark Matter (DM) candidate. This follows from the double suppression of its decay rate, which depends on the 



small R-parity violating couplings and it is suppressed by the Planck scale [10|, LUj □■ Interestingly, gravitino decays 



produce monochromatic photons, opening therefore the possibility to test this scenario with astrophysical searches 



Requiring the model parameters to correctly account for observed neutrino oscillation parameters |13| implies that 
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^ Tlic gravitino as a Warm Dark Matter candidate in the BRpV model with gauge mediation and its collider implications have been 
studied in I12I1 . 
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expected rates for gamma-ray lines produced by gravitino decays of mass above a few GeV would be in conflict with 
the Fermi-LAT satellite observation leading to an upper bound on the gravitino DM mass. The bound on the 

gravitino mass with bilinear R-parity violating couplings holds under the assumption of universality in gaugino masses 
(see [3l and references therein) . Here we study the conditions of non- universality in gaugino masses to relax the the 
gravitino mass constraints. We show how the bound remains even if we assume non-universal gaugino masses, though 
somewhat less stringent. Turning to the implications at collider experiments such as the Large Hadron Collider (LHC) 



these have been discussed in a series of earlier papers |16l - l21| . It is important to stress the expected signatures are 



basically the same already studied in the usual BRpV scenarios. 

In the next Section we briefly introduce the BRpV model, in section lllll we discuss gravitino decays and cosmological 
relic abundance and explain our numerical procedure. In particular we obtain an upper limit for the gravitino mass and 
discuss it both for universal and non- universal gaugino masses. In Section ITVl we briefly discuss collider implications, 
and finally summarize the paper in SecfVl 



II. BILINEAR R-PARITY VIOLATING MODEL 

Here we work in a constrained minimal supergravity model in which the gravitino is the lightest supersymmetric 
particle (LSP). The simplest R-parity violation scenario is assumed, in which the superpotential contains bilinear 
R-parity violating terms 0, Q 

W = WMSSM+e^L^Hu, (1) 

where Wmssm is the superpotential of the minimal supersymmetric standard model (MSSM) and the parameters 
Ci characterize the bilinear R-parity violation, with the flavour index i = 1,2,3 running over the generations. The 
soft supersymmetry breaking Lagrangian contains, in addition to the R-parity conserving operators V^oft^"^^^' ^ term 
associated with the R-parity violation contribution; 

Vsait = KofP^ + BitiLiHu- (2) 

The Ci and Bi terms induce vacuum expectation values Vi for the scalar neutrinos, generating a mass mixing between 
neutrinos and neutralinos. As a result, one finds that neutrino acquires a mass at tree-level given by 0] 

^.ec. ^W+AW^ I^p^ (3) 

where Aq = det(M^p) = —MiM2fJ.'^ + ^iiVdVu{Mig^ + M2g'^) is the determinant of the MSSM neutralino mass matrix 
and Ai = ^Vi + VdCi is the alignment vector. It is worth mentioning that the value of |Ap is almost fixed by the very 
precise determination of the atmospheric scale, mainly by MINOS and other accelerator experiments. The other two 
neutrinos acquire a calculable mass only at the one-loop level Q as required to explain solar and reactor neutrino 
data. The BRpV model provides a scenario where all current neutrino oscillation data can be accounted for, i.e. it can 
accommodate the required values of the neutrino mass squared differences and mixing angles inferred from neutrino 
oscillation studies llSl . We refer the reader to N, M for more details about the model. 



For the rest of the paper, we assume the Constrained MSSM scenario (CMSSM) [2^ where the soft supersymmetry 
breaking parameters mo,mi/2 and Aq are assumed to be universal at the supersymmetric Grand Unification (GUT) 
scale. Thus, the model depends upon the following eleven free parameters: 

mo, mi/2, tan/3, sign(/x), Aq, ei, Ai. (4) 



Here, mi/2 and mo are the common gaugino mass and scalar soft SUSY breaking masses at the unification scale, 
tan /? is the ratio between the Higgs field vacuum expectation values and Aq is the common trilinear term. 
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When supersymmetry is promoted to be a local symmetry of nature, the resulting theory requires a supermultiplct 
which includes the gravitino. After supersymmetry breaking, the gravitino becomes massive via the superhiggs 
mechanism. Depending on details of the underlying supersymmetry breaking mechanism the gravitino mass ttlq 
can lie anywhere between 0{eV) and OiTcV). We take ttiq as a free parameter which does not fix the scale of 
soft-supcrsymmetry breaking parameters 23 1. 



III. GRAVITINO COSMOLOGY 



In this 
Here we 



There have been several studies of gravitino dark matter in R-parity conserving supersymmetry 
case the lightest supersymmetric particle is stable and is potentially a viable dark matter candidate 
consider this issue within the simplest R-parity violating scenario. In order to study the region of parameter space of 
the model which can accomodate neutrino oscillation data, we have performed a numerical analysis using the SPheno 
package 27|, which calculates the renormalization group equations at two loops, generating the full supersymmetric 
particle spectrum. It also includes the one-loop calculations of the neutrino masses in the BRpV model, required in 
order to account for solar neutrino conversion. For a fixed set of CMSSM parameters (mg, TOi/2, tan/3, sign(/x), Aq) 
we determine the set of R-parity breaking parameters and A^ responsible for generating neutrino mass squared 
differences and mixing angles consistent at 3a with the measured values [13|. We repeat this procedure fixing Aq — 
— 100 GeV sign(/.j) = +1 and scanning over the other CMSSM parameters in the range: 



240 < mi/2 < 3000 GeV, 

3 < tan/3 < 50, 
200 < Mo < 1000 GeV. 



(5) 



A. Gravitino dark matter relic density 



Gravitinos are produced in the early Universe after the reheating phase through particles scattering occurring in 
the thermal plasma, the dominant contribution coming from SUSY quantum chromodinamycs processes 28l432l| . The 
gravitino relic abundance critically depends on the reheating temperature Tr and it is given by: 



28 



29 



3m| 



In 



gtiTn) J VlOO GeV/ \10^° GeV 



(6) 



with Mi{Tf;) and gi(Tji) respectively the gaugino mass parameters and gauge coupling constants at Tfj energy scale. 
The index i runs over the Standard Model gauge group factors and the constants uji and ki arc uji = (0.018, 0.044, 0.117) 
and ki = (1.266, 1.312, 1.271). The gaugino masses and the gauge coupling constants can be evaluated at the energy 
scale Tr using the renormalization group equations (RGEs), which at one-loop level are given as 



a,{mz) - — ^In 



-I -1/2 



gi{mz) 



Miimz). 



(7) 
(8) 



In the MSSM, the beta function coefficients are p''^^ — (11,1,-3). Assuming universal gaugino soft masses at the 
supersymmetric GUT scale, their values at the electroweak scale {mz) follow the relations ~ 3.IA/2 — 5.9Mi. 

In figure [T] the black lines show the combination of gravitino masses and reheating temperatures for which the 
gravitino relic abundance is consistent with the dark matter density inferred from astrophysical observations [33[, 
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FIG. 1. Region (white) in the {m^jTii) plane consistent with neutrino oscillations, gravitino dark matter, gamma-ray line 
searches, and gluino searches. In the orange area a consistent gravitino relic abundance would require gluino masses excluded 
by present collider searches. In contrast, the blue area is viable but requires too large gluino masses, Ms > 6000 GeV. The 
yellow region is excluded by astrophysical gamma-ray line searches. 



f^DM^^ = 0.1123. There, we have assumed gaugino universality and the two curves correspond to two different values 
for the soft gluino mass parameter M3. The orange region (upper dark gray) requires values of M3 smaller than 550 
GeV in order to obtain the correct gravitino relic abundance consistent with WMAP. Present bounds on the gluino 
mass from recent searches at the LHC already exclude this region 34" 4^1 • On the other hand, the blue area (lower 
dark grey) is viable but corresponds to M3 > 6000 GeV which, though phenomenologically acceptable, is theoretically 
disfavored if supersymmetry is supposed to "protect" the hierarchical problem. 

We notice that for gravitinos in the mass range up to few GeV, the corresponding values of reheating temperature 
are relatively low and compatible with the lower bounds that can be inferred from cosmic microwave background 
radiation observations [41 1 . 



B. Gravitino decays 



In the presence of R-parity breaking, as in the bilinear R-parity model we have considered above, the LSP decays. In 
particular, if the strength of the bilinear R-parity violating parameters is chosen so as to reproduce the neutrinos masses 
and mixing angles indicated by neutrino oscillation experiments [13| , a neutralino LSP would decay with a lifetime way 
too short when compared with the age of the Universe. However, if the LSP is a gravitino, the double suppression 
provided by the smallness of the R parity violating parameters and the Planck-scale suppression of the coupling 
governing the decay rate greatly increase its lifetime, making it a perfect dark matter candidate lOl lllj. It would 



also provide an example of the generic expectation that gravitational interactions break global symmetries 
in this case R-parity and lepton number. 



42 



4a, 



The unstable gravitino dark matter scenarios can potentially be tested with indirect astrophysical dark matter 
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searches. In the bilinear R-parity breaking model under consideration, the gravitino decays as G — !• 1^7 with a width: 

n3 



r = r(G ^ ^ ya) 



—\u~ 



where the R-parity breaking mixing parameter iC^-yiyp from the 7x7 neutralino mixing matrix, is [l 



a=i+i 



(9) 



(10) 



where the A^-coefficients denote the neutrino projections onto the gauginos. Following 
perturbatively as 



/i^g^ sin^ ( 



'w 



'yv I 



4A2 



(M2-Mi)2|A|2, 



this can be calculated 



(11) 



with Ao = det(A/;^o) = -MiAhfi^ + \^iVdVu{Mig'^ + Ahg''^) and A,; = + Vdti- 

For each set of parameters generated through the scanning procedure given in Eq. ([S]) yielding the correct values 
of the neutrino oscillation parameters, we compute the gravitino lifetime, using equations ^ and (jlip . This decay 
mode is particularly interesting from the point of view of indirect dark matter detection. Indeed monochromatic 
photons of GeV energies are generally not expected to be produced by conventional astrophysical processes. For 
this reason, the detection of gamma-ray lines would be a striking signature of dark matter processes, pointing either 
to annihilations 



44 



or to decays of dark matter particles [51l - l56| . Search of gamma-ray lines have been recently 
performed using the data of the Fermi-LAT satellite |1^ . [isl . The derived upper bounds on the gamma-ray line fluxes 
can be used to constrain the unstable gravitino dark matter model under consideration. 

In figure [5] we present the lower bounds on the gravitino lifetime for dark matter gravitinos decaying into iry 



m 



These constraints have been computed assuming a Navarro-Frenk- White (NFW) dark matter density profile 57[ and 
for the region of observation dubbed as "Halo" in Ref. 15|. The bounds are not too sensitive to the exact shape 
of the dark matter profiles or region of observation considered. At energies below ~ 1 GeV we consider the bounds 
on gamma-ray lines obtained in Ref. 58[ by analyzing the data from EGRET. We traslate the upper limits on the 
gamma-ray line fluxes into bounds on the gravitino lifetime for gravitino decays into 1^7 [58[. We consider a NFW 
density distribution while a shallower isothermal profile would lead to a bound a factor two less stringent. Finally, 
we note that gravitino decays into three body final states could be relevant for gravitino masses larger than those 
required in our case 
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The area between the two black lines in figure [5] corresponds to the region of the parameters compatible with 
neutrino physics. We notice that the two curves correspond approximately to constant values of 77ii/2 = 240 GeV 
(lower line) and mi 12 = 3000 GeV (upper line). Indeed, once the constraints from neutrino oscillations are imposed, 
the matrix element iCZ-Yjyp determining the gravitino lifetime, depends mostly on mi/2- Universal gaugino masses of 
TOi/2 = 240 GeV and mi/2 = 3000 GeV lead at the scale Mz to M3 k. 550 GeV and M3 « 6000 GeV respectively. 

Taking into account the Fermi-LAT and EGRET bounds on gamma-ray lines from dark matter decay (yellow region 
in figure [5]) and assuming the gravitino as dark matter particle, we can derive an upper bound on the gravitino mass 
of the order ^ 2 GeV. Thus, the grey area in figure [5] corresponds to the region of the parameter space which 
simultaneusly explains neutrino oscillation data and satisfies the constraints from gamma-ray line searches. 

We now translate the bounds from neutrino oscillation physics and gamma-ray line searches to the (m^^Tji) plane 
shown in figure [TJ They correspond to the yellow area (light grey). We see that assuming the gravitino as dark 
matter candidate in BRpV and imposing the constraints from neutrino oscillations and gamma-ray line searches we 
can derive an upper bound on reheating temperature of the order Tj^ ~ 10® GeV and an upper bound on the gravitino 



mass of the order 



2 GeV. Similar results are expected in other R-parity violation schemes, such as considered 



in Ref. 



55|. 
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FIG. 2. Allowed gravitino mass-lifetime region (grey color) consistent with neutrino oscillation data and astrophysical bounds 
on gamma-ray lines from dark matter decay. The yellow region is excluded by gamma-ray line searches (Fermi and EGRET 
constraints are respectively above and below 1 GeV). The lower and upper black lines correspond to mi/2 ~ 240 and 3000 GeV 
respectively. 



C. Non- universal gaugino masses 

We now study the effects of non-universal gaugino masses on the gravitino mass upper bound. The non-universaUty 
effects enter in the gravitino lifetime through the neutrino-^hotino" mixing parameter |J7^i/p. At the unification 
scale non-universal gaugino masses can be parametrized as j6ll | 

Ma = 171,^^(1 +Sa), (12) 

with the parameters 6a, a ~ 1,2,3, characterizing the deviation from universality. For illustratation we choose the 
ranges Si^2 = (—1, 1), keeping 63 = 0, and fix a typical CMSSM point satisfying all the phenomenological constraints 
with mi/2 = 500 GcV, niQ = 1000 GeV, Aq = -100, tan/3 = 10, and sgn{fi) > 0. Then we use SPhcno with a 
random set of da values to calculate |J7^i/p, with the best possible fit to neutrino masses and mixings for each point 
of the scan. We also check that limits on sparticle searches are obeyed, e. g. m-± > 103 GeV, and that a non-bino 



neutralino has a mass larger than 50 GeV [6l|, |62| . The results are shown in the yellow (light gray) region of figure [3] 
where the ratio M2/M1 is calculated at the clcctrowcak scale. In the green (dark gray) region the mass neutralino is 
larger than 50 GeV, while in the solid black line 62 = 0. 

From the approximate expression for jC/^i/p in cq. PT|) one would expect vanishing values at M2 » Mi. However, 
from the loop-corrected neutralino mass matrix calculated from SPheno, one obtains that |C^7i/p has a minimum 
non-zero value at M2 ~ Mi . In figure [3l the dashed lines indicate the point where the gaugino masses arises from 
universal conditions. From the gravitino lifetime in (jH]), one sees that for neutralino masses larger (smaller) than the 
CMSSM reference value, the gravitino mass bound is weaker (stronger). For example, for mi/2 ~ 500 GeV, we have 
|J7^^P « 1.6 X 10"". From cq. ^ the maximum gravitino mass can be expressed as 
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M„{Agut) =mi/2(l + <S<.), 




10' 

MJM, 



FIG. 3. Neutrino-photino mixing |f/-y^|^ as a function of the low energy ratio M-z/Mi when Ma ~ rn-i/2{^ + 5a) have been 
assumed at the unification scale (with So, — 0)). We have set = 500 GeV, mo = 1000 GeV, Aq = —100 GeV, tan /? = 10 

and sign(/i) — +1. In the green (dark gray) region m^o > 50 GeV. The explicit minimum of \Ujv\^ is shown in the zoomed 
gray area. 



where r^*" is the minimum gravitiiio lifetime allowed by gamma-ray line searches (see figure [2]). When M2 ~ Mi, 
we obtain the minimum value for iC/^zyp ~ 4 x 10"^'', and an upper bound on gravitino mass of order itIq ^ 7 GeV. 
This illustrates the relative robustness of the gravitino mass bound against deviations from gaugino universality. 

While this holds for a given CMSSM point, we have not been able to find other points in parameter space where 
the gravitino mass bound changes by more than an order of magnitude. The point is that, even though radiative 
effects in the neutralino mass matrix may change by three orders of magnitude (see figure [3]) the bound changes only 
as the 1/3 power of that, according to eq. (|T3l) . hence is relative stable. 

Let's now briefly comment about implications on Big Bang Nucleosynthesis (BBN). In R-parity conserving scenarios 
with gravitino dark matter, the neutralino has a large lifetime since its decays into the LSP are suppressed by the 
Planck scale thus it may decay during the Big Bang Nucleosynthesys epoch, spoiling its predictions [63| (BBN demands 
a NLSP lifetime less than 0.1 s [fi^). In contrast, in the model under consideration, the next to lightest supersymmetric 
particle (NLSP) decays occur well before the BBN epoch because of the presence of the gravity unsuppresscd R-parity 
violating interactions, keeping therefore the successful BBN predictions of the light element abundances. 



IV. PROSPECTS FOR COLLIDER SEACHES 



When R-parity is conserved, all supersymmetric particles undergo cascade decays to the next to lightest supersym- 
metric particle, which subsequently decays (of course, with gravitational strength) to the gravitino. The implications 
for coUider searches and cosmology strongly depend on which superpartner is the NLSP. In BRpV models, in addition 
to generating the neutrino masses, the neutralino-neutrino mixing also induces NLSP decays into Standard Model 
particles, strongly correlated with the neutrino oscillation parameters 16-18[. Since R-parity violating couplings 
are not so small, displaced vertices are expected in the NLSP decay 19l421|. In what foUows we will consider the 
neutrahno as the NLSP. 
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In R-parity conserving scenarios the neutralino as the NLSP has the following decay channels 64 1 



Xi ^ ZG, 



(14) 



In the presence of R-parity breaking additional decay channels exist 18[ , namely, 



x? 



(15) 



Neutralino can also decay to three fermions by scalar quark and scalar lepton exchange in R-parity violating models. 
The three channels in (|14p are Planck-mass-suppressed and, for the gravitino mass range of interest, are negligible 

Indeed, Br{xi 



-fG) < 10^5 for > 10 keV 



12[. The decay into the Higgs 



compared with those in cq. (|15p 
boson is scalar mixing suppressed, while the radiative channel is loop suppressed. As a result, decays to gauge bosons 
are dominant for large mp. Therefore these collider signals are independent of the fact that the gravitino is lightest 
supersymmetric particle or not. In particular, the predictions at colliders for a neutralino LSP in the CMSSM with 
BRpV studied in 16l421j. such as the displaced vertex signals illustrated in figure SI remain unchanged. 





0.001 



400 



m^o{GeV) 



FIG. 4. Neutralino NLSP decay length as a function of its mass. For illustration we show the results of a scan with Aq 
GeV, sign(^) = -1-1, tan/3 = 10, 200 < Mo < 1000 GeV and 240 < mj/a < 1000 GeV. 
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V. SUMMARY 

We have considered the CMSSM model with bilinear R-parity violation with gravitino as LSP. By imposing the 
constraints from neutrino oscillation data, dark matter relic density and gamma-ray line searches at Fermi-LAT and 
EGRET we have shown that the gravitino does provide a viable radiatively decaying dark matter particle, provided 
its mass and reheat temperature are bounded as < 1-10 GeV and < few xlO'' — 10^ GeV (as we saw the 
bounds get looser if the universality hypothesis in gaugino masses is relaxed). The expected signatures associated 
to the NLSP decays at collider experiments like the Large Hadron Collider do not depend on the presence of the 
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gravitino and so arc the same as those previously studied 16l-l2ll| . In particular neutrino mass and mixing parameters 
may be reconstructed at accelerator experiments by measuring the ratio of semileptonic neutralino decays branching 
ratios induced by the charged current. 
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